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Abstract 

Polypyrrole is an attractive polymer for use as a high-energy-density secondary 
battery because of its potential as an inexpensive, lightweight, and noncorrosive 
electrode material. A mathematical model to simulate cyclic voltammograms for 
polypyrrole is presented here. The model is for a conductive porous electrode film on a 
rotating disk electrode (RDE) and is used to predict the spatial and time dependence of 
concentration, overpotential, and stored charge profiles within a polypyrrole film. The 
model includes both faradaic and capacitive charge components in the total current 
density expression. 


Introduction 

Electrochemically synthesized conducting polypyrrole has received much attention 
as a secondary battery electrode because it has high conductivity, high specific charge 
capacity, and stable and reversible electrochemical redox behavior (1-4). Polypyrrole, 
in its conducting polycationic form (PPy + ), is readily obtained by a simple one- 
step electrochemical oxidation of the pyrrole monomer. This reaction is electronically 
irreversible and fairly complex. A detailed discussion of the mechanism of synthesis is 
beyond the scope of this study. 

Electronically conductive polypyrrole has a number of potential technological 
applications. It is used 1) as a battery electrode, 2) in display devices, 3) in electronic 
devices, and 4) as a component of new high-energy-density secondary batteries. All 
of these applications rely on the ability of the polypyrrole to switch between the 
conductive and nonconductive states. Polypyrrole in its conducting polycationic form 
can be transformed by cathodic reduction into the neutral insulating state (PPy°) with 
conductivity less than 10~ 5 ohm -1 cm -1 . This film is extremely unstable and easily 
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reoxidized to yield PPy + . In contrast, oxidized polypyrrole film is stable in a wide 
range of solvents. 

An experimental cyclic voltammogram of polypyrrole approximates the behavior 
of a quasi-reversible couple with the distinct characteristic of a large capacitive 
background current (5). Some authors (6-8) have tried to analyze the capacitance 
effect in order to separate the faradaic current from the capacitive current. In spite 
of this, the capacitance effect is still not fully understood. Feldberg (6) suggested 
a mathematical model of the capacitive current with the electrochemical switching 
polymer films, but failed to include transport effects (diffusion and migration) of the 
counter-ions. Pickup and Osteryoung (7) developed a mathematical model for potential 
step chronoamperometry and concluded that a polymer film could be modeled as a 


porous electrode. The porous electrode model they used is not satisfactory because it 
is based on the assumption that faradaic reactions are negligible. 

The model presented here is used to simulate potentiostatic cyclic voltammograms 
for a polypyrrole film on a RDE. The rotating disk electrode was chosen because of 
its well-defined hydrodynamics. The model consists of a combination of the models 
presented by Feldberg (6) for the capacitance effect on the current density and by Ryan 
ei al. (9) for a conductive porous layer on a rotating disk electrode. The polypyrrole 
film is treated as a porous electrode with a high surface-to- volume ratio and a large 
double layer capacitance, which is proportional to the amount of oxidized polypyrrole 
film (4-7). These model conditions were chosen so that the simulated predictions could 
be compared to experimental results available in the literature. The faradaic and 


capacitive current components associated with the electrochemical switching of a film 
between the insulating and electronically conducting states are included in the model. 

The transport equations of the model are based on the conservation of mass and 
charged species. The equations include migration of charged species in an electric field, 
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diffusion of charged and uncharged species, and the electrochemical reaction that occurs 
within the porous polypyrrole film. Since the surface of a solid electrode is composed 
of active (oxidized) and inactive (reduced) sites during charge and discharge processes, 
averaged values are used to describe local variables within a volume element throughout 
the porous layer to account for the effect of nonhomogeneity of the electrode surface 
(9). 

The model is used to predict concentration, overpotential, and stored charge 
profiles within the polypyrrole film as a function of position and time. The model 
is also used to predict the dependence of the faradaic, capacitive, and total current 
density on applied overpotential. 


Model Development 


The model presented here is for predicting potentiostatically controlled cyclic 
voltammograms for a porous polypyrrole film on a RDE in a one-compartment cell. 
The one-compartment electrochemical cell contains a platinum RDE coated with 
polypyrrole film, a platinum counter electrode, and a saturated calomel reference 
electrode (SCE). A schematic diagram of a typical experimental system is shown in 
Fig. 1. The Luggin capillary tip of the reference electrode is placed as close as possible 
to the center of the disk. This enables one to use the reference electrode to detect the 
solution potential near the working electrode. The electrolyte consists of 0.1 M LiC 104 
in acetonitrile. In the solution, LiC 104 dissociates into charged Li + and C10 4 . 

During a potential sweep, a polypyrrole film can be cycled repetitively between its 
oxidized and reduced forms in the absence of oxygen with no evidence of polypyrrole 
decomposition (4). The reaction occurring during the potential sweep is assumed to 
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be described as follows: 

PPy° + C10 4 - ^ PPy+C10 4 - + e“ [1] 

The rate of this reaction is controlled primarily by the available electroactive area and 
the transport rate (diffusion and migration) of CIO^. 

Figure 2 is a schematic representation of the regions close to the RDE surface that 
are relevant to the development of the model equations. As shown in Fig. 2, there 
are two separate regions: the porous electrode region of width £ PE and the solution 
diffusion region of width S RE . 

The dependent variables in the diffusion layer are: the concentration of Li + (c+ ), 
the concentration of CIO^ (c_), and the potential of the solution phase ($ 2 )- In the 
porous polypyrrole region, the dependent variables include the ones in the diffusion 
layer, the local faradaic charge per unit volume ( Q F ), and the potential of the solid 
polymer phase ($ 1 ). Values for these unknowns depend on the perpendicular distance 
from the electrode surface ( y ) and time (f), and they are obtained by solving the 
governing equations with associated boundary conditions. 

Assumptions 

The assumptions used in developing the governing equations for the model are as 
follows: 

1) no homogeneous reactions occur in the bulk solution. 

2) the solution is Newtonian and isothermal. 

3) transport occurs only in the axial direction (i.e., y direction). 

4) the platinum current collector is solid and inert. 

5) the differential capacitance of the platinum is negligible. 

6) the differential capacitance of the double layer is proportional to the amount of 

oxidized polypyrrole. 
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7) the differential capacitance of polypyrrole is independent of potential. 

8) the electrochemical reaction occuring within the porous electrode is a pseudo- 
homogeneous reaction whose rate is based on the Butler- Volmer equation, as 
presented below. 


Governing Equations — Solution Diffusion Layer 

Mass transfer in the diffusion layer is governed by the following material balance 
equation for species i: 

dc ; 

— = -V . N ; where i = +,- [2] 

Eq. [2] depends on the flux, N;, of species i which is due to migration in the electric 
field, diffusion in the concentration gradient, and convection in the flow field 


Nj = — z\U\ i r ’cjV $2 — AVci + CjV where i = +,— 


[3] 


The ionic mobility, iq, is assumed to be described by the Nernst-Einstein equation (10) 

D\ 


U\ = 


RT 


where i = +, 


[4] 


Only the axial component of Eq. [3] is considered in the model and the velocity 
component in that direction depends only on the normal distance from the electrode 
surface according to the no-slip condition (9) 

v y = -a'Sl^(y - y PE ) 2 [5] 

Combining Eq. [2], [3], [4], and [5] yields the governing equation for species i within the 
solution diffusion layer on the RDE: 

da Z;D;F( , dc; 0$ 2 \ n d 2 C, [fi . , 

m = n^{ c> -w + dy-^) + V where , = + ’- |6] 

The electroneutrality condition provides the additional governing equation needed to 
solve for $2 within the diffusion layer 

y: Z\C\ = z + c+ + Z—C— = 0 [7] 



7 


Governing Equations — Porous Electrode Layer 

To account for the nonhomogeneous structure of the polymer film, macroscopic 
properties are used to describe the porous electrode layer in terms of measurable 
parameters without regard for the actual geometric detail of the pore structure. Two of 
these properties are the porosity (e) and the MacMullin number (JV M PB ). The porosity 
represents the void volume occupied by the electrolyte per unit volume of porous region. 
The MacMullin number (11) is defined as the ratio of the tortuosity to the porosity: 

A m,pe = T / c [8] 

The porosity and MacMullin number are assumed to be constant. 

Within the porous electrode layer, the material balance equation is formulated in 
terms of average quantities that include the reaction rate term for the electrochemical 
reaction 

&€C' 

= -V • N; + R[ where i = +, - [9] 

oi 

and R[ is the production rate of species i due to a pseudo-homogeneous reaction 
(electrochemical reaction) within the porous electrode. 

The flux expression, Eq. [3], is modified to account for the porous and tortuous 
structure of the polypyrrole film. This is accomplished by replacing the diffusion 
coefficient in Eq. [3] by an effective diffusion coefficient within the porous layer 

A,e = ~77 ~ ~ where i = +,- [10] 

^M,PB 

Thus, the flux of species i within the solution phase of the porous layer becomes 

N; = — Z|I?1 p CiV $2 ~ Vq where i = +, - [11] 

iV M,PE 

and where it has been assumed that there is no convective flow within the porous 
layer. Combining Eq. [11] with the material balance expression, Eq. [9], subject to the 
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assumptions mentioned above, yields the governing equation for species i within the 
porous electrode layer: 

chcj ZjDjF ( , 6*9*^ A &Ci . _ ri 

dt N U „RT\ C ‘ dy* + dy dy ) + JV„.„ dy* + 1 W ‘ ’ [ 1 

The production rate of species i due to an electrochemical reaction is given by 

( 12 ): 

i?| = -^p F where i = +, - [13] 

and j p is the faradaic transfer current per unit volume (A/cm 3 ) of the porous region 
and is disscussed later. 

It is a consequence of the assumption of electroneutrality that the charge which 
leaves the solid phase must enter the solution phase. This is expressed as (12) 


V-t 1 +V.t 2 = 0 


[14] 


where ii is the superficial current density (current per projected area of the electrode) 
in the solid electrode phase, and *2 is the superficial current density in the solution 
phase. In the solid phase, the movement of electrons is governed by Ohm’s law 


ii = ~<r(l - e) 


dy 


[15] 


wVi*»r#=> rr ic rondiirtivit.v nf nnlvnvrrnlp anr) sr(1 — pi is f>nni valpnt to f»ffprt.ivf» 


r — j rj 


solid conductivity. It is assumed that the conductivity of polypyrrole is proportional 
to the degree of the oxidation of polypyrrole as follows: 

[16] 


^MIN ^MAX 


Q* 


Q* 


(Qf - Q F,M AX ) + °MAX 


'F,MIN v F,M AX 

where <r MIN represents the conductivity of the polypyrrole film at minimal charge state 
(Q F MIN ), and <r MAX represents the conductivity of the film at its maximal charge state 

( O Vi /•iii'rnnf on ~\r itj + Vi a caIii! mn T"\ Vi sj c #=» ic fTi Vitr 

V^F MAX/* J- lie Linivnt HIV JVlUVlVll • vil 


*2 = 


N, 


K d$2 
m,pe dy 


[17] 
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where k. is the specific conductivity of electrolyte and is assumed to be constant. The 
term k/N m pe is the effective solution conductivity. Substitution of Eq. [15] and [17] 
into the conservation of charge equation, Eq. [14], yields 


V • 


-<r(l - e) 


ggi 

dy 


+ V 


K d$2 \ 
n m,pb &y ) 


= o 


[18] 


Electroneutrality, Eq. [7], is also used as the governing equation for $2 in the porous 


region. 


The rate of accumulation of the faradaic charge contained in the oxidized polypyr- 


role per unit volume, Q F , is assumed to be related to the faradaic transfer current, j F , 
as follows: 

dQ F _ . 

di ~ Jf 



Current Density 


The faradaic transfer current per unit volume, j F , is assumed to be given by the 
Butler- Volmer equation (13): 



[ 20 ] 

where a is the electroactive surface area per unit volume, i Q REF is the exchange current 
density at a given reference concentration (c i REF ), and r/ is the overpotential. The 
overpotential is defined as: 


7 /=($i -$ rb )-(*2 [ 21 ] 

where 4> RE is the solution potential at the Luggin tip, and f/ REF is the open circuit 
potential at a given reference concentration (c ; REP ) and is relative to SCE. Anodic and 
cathodic current densities are taken to be positive and negative, respectively. Note 
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also that o: a + a c = n. It can be seen that the local transfer current predicted by 
the Butler- Volmer kinetic expression, Eq. [20], depends upon the difference between 
the potential of the solid phase and that of the adjacent solution within the porous 
electrode. 

The anodic faradaic current in the porous polymer film leads to charging of the 
double layer within the pores of the polymer film in a manner consistent with that 
proposed by Feldberg (6). That is, the amount of capacitive charge that goes to 
charging the double layers within the pores of the porous film, Q c , is related to the 
amount of the faradaic charge added to the polymer film by the faradaic reaction: 

Q c = a(ri - rj PZC )Q F [22] 

where a* is a proportional constant which is assumed to be independent of potential, 
and r] PZC is the total overpotential across the double layer at the point of zero charge 
(PZC) which is given by 

W = { (*1 - *RE ) - (*2 - *RB ) - t'RBF } [23] 

t ' PZC 

The capacitive transfer current per unit volume, j c , associated with charging of 
the double layer is defined as follows: 

. _dQ c 
Jc dt 

Substituting Eq. [19] and [22] into Eq. [24] shows that 

j c = a* (q f ^ + {v- Vpzc)j^j [25] 

with the assumption that $ RE = ($ RE ) PZC . The total transfer current per unit volume, 
j T , is defined as the sum of the faradaic and capacitive transfer currents per unit 
volume: 

Jt 3y "b 3c 



[26] 
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Eq. [26] can be used to predict the total current density, i T (current per projected 
area), by integrating j T over the porous layer: 

rv=y pe 

* T = / JT d y [27] 

J y=0 

Note that the faradaic current density, i F , and the capacitive current density, z c , can 
be obtained in a similar manner from j p and j c , respectively. 

Boundary Conditions 


The boundary conditions at the polymer film/substrate interface ( y = 0) are 


ziD>F d$ 2 

- RT Ci dy - Di dy = ° = 

[28] 

— E + $ 

^APP ' * RE 

[29] 

ZjCj = Z+C+ + z_c_ = 0 
i 

[30] 

where E APP and $ RE are constants. The boundary conditions at 
film/diffusion layer interface ( y = y PE ) are 

the polymer 

N^y [porous layer = Njj [diffusion layer where i = T , 

T = ° 

dy 

[31] 

[32] 

ZiC, = z+c+ + z-C- - 0 
i 

The boundary conditions in the bulk solution (y = y RE ) are 

[33] 

Ci = c iREF where i = +,- 

[34] 

$2 = $ RE 

[35] 
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Initial Conditions 

The initial concentration of each species i is set equal to its bulk concentration: 

c i = c i,REP where i = +,- [36] 

For convenience, it is assumed that the polypyrrole film is initially fully reduced and is 
ready to be oxidized. Consequently, the faradaic charge per unit volume in the porous 
region is initially equal to Q F MIN , a minimal charge state: 

Qf = Qp.min [37] 


Solution Method 

The governing equations and boundary conditions can be written in finite difference 
form and solved by using Newman’s pentadiagonal block matrix equation solver (14). 
Implicit stepping can be used for the time derivatives (15). In the diffusion layer, 
Q F and $2 are dummy variables treated as constants and are set arbitrarily equal to 
zero. During the potential sweep, the applied potential (£ APP ), which is the difference 
between $1 at y = 0 and $2 at y=y RE (^re^)’ * s changed linearly between E im and 
-Z?fin according to a specified scan rate (u s ). 

Results and Discussion 

Cyclic voltammograms of polypyrrole can be predicted using the model described 
above. The values used for the fixed parameter are given in Table I. Figure 3 
shows predicted cyclic voltammograms for the polypyrrole film at various scan rates. 
Comparison of these predictions with those obtained experimentally for a similar 
system, as shown in Fig. 4, shows qualitative agreement. 
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Figure 5 shows the total current density and its components at a 100 mV/sec 
scan rate. Note that when the maximal quantity of oxidized polypyrrole is formed, no 
further oxidation of the film occurs and the current is entirely dominated by capacitive 
charge. 

Figure 6 shows the dynamic profile of the concentration of CIO^ in the porous 
electrode region at an 100 mV/sec scan rate. For convenience, the concentration of 
CIO^ was made dimensionless relative to its reference concentrations (c_ REF ). The 
position coordinate was made dimensionless by using the thickness of the film (6 PE ) 
so that £ = 1 represents the interface between the diffusion layer and the polymer 
film. Time was made dimensionless by using the time required to scan the applied 
potential from E mi to £ FIN (r°) so that £ = 1 represents the intermediate point where 
the polypyrrole film is totally oxidized. 

During the oxidation process, the reacting species (C10^~) is transported from the 
bulk to the porous electrode layer where it diffuses and migrates to reactive sites within 
the porous layer. During reduction, the opposite is true. Since the effective diffusivities 
of Li + and C10^~ within the porous layer are smaller than the free stream diffusivity 
of these species, then concentration gradients within the porous region must be larger 
to make up for the slower movement of the ions. The concentration profiles of both 


species are the same because of the electroneutrality condition. 

Figure 7 shows the faradaic charge per unit volume due to the electrochemical 
reaction within the film at a sweep rat* of 100 mV/sec. The faradaic charge per unit 
volume was made dimensionless by using the maximum faradaic charge value (Q FiMAX ) 
as the reference point. The rate of charge accumulation is faster in the outer layer of 
the film during oxidation because of diffusion and migration effects. Note that charge 
accumulates rapidly throughout the entire film. 
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Summary 

Theoretical calculations were carried out that describe quantitatively the current 
responses of the quasi-reversible behavior of electronically conductive polypyrrole. 
Porous electrode theory was applied to the conductive polymer to provide the basis for 
understanding the transport behavior of the switching process within the polymer film. 
The effects of diffusion and migration were included because it is likely that counter- 
ion movement limits the oxidation and reduction rates. The net cyclic voltammetric 
current was decomposed into two components: a capacitive current due to the double 
layer charging and the faradaic current due to an electrochemical reaction. It may be 
possible to use this model together with experimental data and parameter estimation 
techniques to determine the transport and kinetic parameters of the model. 
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NOTATION 

a specific surface area of the porous material, cm -1 

a 0.51023 

a* double layer constant, V -1 

c+ concentration of Li + , mol/cm 3 

c_ concentration of CIO^, mol/cm 3 

Cj concentration of species i, mol/cm 3 

cf dimensionless concentration of species i, (Cj/c iRBF ) 

c i,ref initial concentration of species i in the porous region, mol/cm 3 

D\ diffusion coefficient of species i, cm 2 /sec 

D\ e effective diffusion coefficient of species i in a porous layer, cm 2 /sec 
E APP applied potential (potential difference between the current collector and 

*»), v 

E Fm final applied potential, V 

E im initial applied potential, V 

F Faraday’s constant, 96487 C/mol 

i c capacitive current density, A /cm 2 

i F faradaic current density, A/cm 2 

i 0 REF exchange current density at reference concentrations, A/cm 2 
i T total current density, A/cm 2 

i\ superficial current density in the porous electrode phase, A/cm 2 

*2 superficial current density in the solution phase, A/cm 2 

j c capacitive transfer current per unit volume, A/cm 3 

j F faradaic transfer current per unit volume, A/cm 3 

j T total transfer current per unit volume, A/cm 3 
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NOTATION (Continued) 


n number of electrons transferred 

Nj flux vector of species i, mol/cm 2 -sec 

Nj jY y component of flux vector of species i, mol/cm 2 -sec 

7V M PE MacMullin number for the porous electrode layer 

Q c capacitive charge of double layer per unit volume, C/cm 3 

Q f faradaic charge of polymer film per unit volume, C/cm 3 

<2* dimensionless faradaic charge of polymer film, (Q F /Q P|MAX ) 

Q F max maximum faradaic charge of polymer film per unit volume, C/cm 3 
Q F MIN minimum faradaic charge of polymer film per unit volume, C/cm 3 
R universal gas constant, 8.3143 J/mol-K 

i?| pseudohomogenous reaction rate of species i, mol/cm 3 -sec 

s\ stoichiometric coefficient of species i 

t time, sec 

T absolute temperature, K 

u\ mobility of species i, mol-cm 2 /J-sec 


U KEF 

V 


y 

VpE 

Vre 

Z'l 


theoretical open-circuit potential at reference concentration, V 


electrolyte velocity vector, cm/sec 
scan rate, V/sec 

electrolyte velocity in the normal direction, cm/sec 

perpendicular distance from a porous electrode/substrate interface, cm 
position of a porous region/ diffusion layer interface in y coordinate, cm 
position of a bulk solution in y coordinate, cm 
proton charge number of species i 



NOTATION (Continued) 


Greek Symbols 

a a anodic transfer coefficient 

a c cathodic transfer coefficient 

6 pe thickness of polypyrrole film, cm 

8 rb thickness of diffusion layer, cm 

e porosity or void volume fraction 

( dimensionless time, (t/r°) 

rj overpotential, V 

Vpzc overpotential at point of zero charge, V 

k solution conductivity, fl -1 cm -1 

v kinematic viscosity, cm 2 /sec 

£ dimensionless distance, (1 — y/S PB ) 

p 0 pure solvent density, kg/ cm 3 

a conductivity of polypyrrole film, fl _1 cm _1 

<r M ax conductivity of totally oxidized polypyrrole film, fl -1 cm“ 

^min conductivity of totally reduced polypyrrole film, fl _1 cm _] 

t tortuosity of porous material 

t° time required to scan the potential from 2£ INI to £’ PIN , sec 

potential in solid phase, V 
$2 potential in solution phase, V 

$ RE potential in the solution at y = y RE , V 
disk rotation velocity, rad/sec 
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Table I. Fixed parameter values used for cyclic voltammograms 

Operating Variables 

T 

298.15 K 


n 

377 rad/sec 


2/re 

0.01 cm a 


^RE 

0.0 V 


E INI 

-0.8 V (vs. SCE) b 



0.4 V (vs. SCE) b 


Electrolyte Properties 

D+ 

1.03 x 10 cmr/se< 


£>_ 

1.81 x 10 -5 cm 2 /se( 


jREP 

1.547 x 10~ 4 mol/ci 


C — jREF 

1.547 x 10~ 4 mol/ci 


K 

0.0165 c 


V 

6.6 x 10 -3 cm 2 /sec 


P 

0.777 x 10 -3 kg/cm 


Electrode Properties 


20 nm 


e 

0.1 


N 

1 * M, PE 

1.0 x 10 6 


Q F.MAX 

344.65 C/cm 3 b 


Qf.MIN 

l.OxlO -14 C/cm 3 


^MAX 

100.0 n-'cm- 1 b 



l.OxlO -14 Q -1 cm 
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Fig. 1. Schematic of a one-compartment electrochemical cell with RDE. 

Fig. 2. Schematic of the surface and associated diffusion layer of the RDE. 

Fig. 3. Simulated current density profiles at various scan rates. 

Fig. 4. Experimental cyclic voltammograms of a similar system (Ref. (4)). 

Fig. 5. Simulated total current density and its components at v s = 100 mV/sec. 
Fig. 6. Dimensionless concentration profile of C10^~. 

Fig. 7. Dimensionless charge profile of polypyrrole film. 
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